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Transverse momentum spectra of charged hadrons with pr < 8 GeV/c and neutral pions with 
Pt < 10 GeV/c have been measured at midrapidity by the PHENIX experiment at RHIC in d+Au 
collisions at ^/s NN = 200 GeV. The measured yields are compared to those in p+p collisions at the 
same ^s NN scaled up by the number of underlying nucleon-nucleon collisions in d+Au. The yield 
ratio does not show the suppression observed in central Au+Au collisions at RHIC. Instead, there 
is a small enhancement in the yield of high momentum particles. 

PACS numbers: 25.75.Dw 



High transverse momentum (px > 2 GeV/c) hadrons They arise from fragmentation of quarks and gluons (par- 
provide an excellent probe of the high energy density tons) scattered with large momentum transfer, Q 2 , in the 
matter created in relativistic heavy ion collisions 1- initial parton-parton interactions 3]. In the absence of 



3 



medium effects, these hard scattering yields in nucleus- 
nucleus collisions should scale with the average number of 
inelastic nuclcon-nucleon collisions N co u (binary scaling) . 
One of the most intriguing observations from experiments 
at the Relativistic Heavy Ion Collider (RHIC) is the large 
suppression of high pr neutral pion and charged hadron 
yields in central Au+Au collisions with respect to p+p 
results scaled by the number of binary nucleon-nucleon 

collisions!! H 00. 

Theoretical studies of parton propagation in high den- 
sity matter suggest that partons lose a significant fraction 
of their energy through gluon bremsstrahlung El Q • re- 
ducing the parton momentum and depleting the yield 
of high p T hadrons H El El El E! This is a 
final-state effect in the spatially extended medium cre- 
ated in A+A collisions. Initial state effects include nu- 
clear modifications to the parton momentum distribu- 
tions (structure functions), and soft scatterings of the in- 
coming parton prior to its hard scattering. These should 
be present in p+A, d+A and A+A. Interpretations of 
Au+Au collisions based on initial-state parton satura- 
tion effects E3 or final-state hadronic interactions E3 
also predict a considerable suppression of the hadron pro- 
duction at high pt- It is therefore of paramount interest 
to determine experimentally the modification, if any, of 
high pt hadron yields due to initial state nuclear effects 
for a system in which a hot, dense medium is not pro- 
duced in the final state. This Letter reports on charged 
hadron and ir° production at midrapidity obtained by the 
PHENIX [il| experiment at RHIC in d+Au collisions at 
ys^" = 200 GeV. The results are compared to those in 
p+p and Au+Au 0, 0, 0] at the same nucleon-nucleon 
center-of-mass e nergy. Similar d+Au measurements are 
also reported in [LjJ, |20j . 

The data were collected under two different trigger con- 
ditions. Minimum bias events with vertex position along 
the beam axis within \z\ < 30 cm were triggered by the 
Beam-Beam Counters (BBC) Ell which cover |^| = 3.0- 
3.9. The minimum bias trigger accepts (88 ± 4%) of 
all inelastic d+Au collisions that satisfy the vertex con- 
dition. A total of 1.2 x 10 7 and 1.7 x 10 7 events were 
analyzed for charged hadron and ir° spectra, respectively. 
A second "photon triggered" sample, requiring showers 
above 2.5 GeV for lead-scintillator (PbSc) and 3.5 GeV 
for le ad-g lass (PbGl) Electromagnetic Calorimeters (EM- 
Cal) [16j in addition to the BBC requirement, is used to 
extend the 7r° measurements to higher px- This trigger 
sampled a total of 1.7 x 10 9 events. 

Neutral pions were measured by the PHENIX EMCal 
via the ir° — » 77 decay. The EMCal covers \r}\ < 0.35 
in pseudorapidity and A(f> = n in azimuth. It con- 
sists of 6 PbSc and 2 PbGl sectors, each covering 77/8 
in azimuth. The data from PbSc and PbGl were ana- 
lyzed separately. The energy calibration is obtained from 
beam tests, cosmic rays, minimum ionizing energy peaks 
of charged hadrons, and the energy/momentum ratio of 



TABLE I: Systematic errors in percent on n invariant yields 
for PbSc (PbGl), as a function of pr (in GeV/c). There are 
3 categories of uncertainty: Type A is a point-to-point er- 
ror uncorrelated between pr bins, type B is pr correlated, 
all points move in the same direction but not by the same 
factor, while in type C all points move by the same factor 
independent of pr- 





type 


PT=2 


Pt=6 


p T =10 


peak extraction 


A 


5.0(5.0) 


5.0(5.0) 


5.0(5.0) 


geom. accept. 


B 


3.0(3.0) 


2.0(2.0) 


2.0(2.0) 


7T° reconstr.eff. 


B 


4.0(4.0) 


4.0(4.0) 


4.5(4.5) 


energy scale 


B 


4.0(4.0) 


9.0(9.0) 


11.0(11.0) 


trigger eff. 


B 




5.0(10.0) 


3.0(3.0) 


trigger norm. 


C 




5.0(5.0) 


5.0(5.0) 


conversion corr. 


c 


2.8(2.8) 


2.8(2.8) 


2.8(2.8) 


total error 




8.6(8.6) 


14(16) 


15(15) 



electrons identified with the Ring Imaging Cerenkov De- 
tector (RICH). The combined uncertainty on the energy 
scale and linearity is <1.5%, determined from the re- 
sponse to identified electrons, and confirmed by the po- 
sitions and widths of the observed 7T° mass peaks. 

Photon-like energy clusters in the calorimeter are se- 
lected via shower profile cuts. The invariant mass for all 
photon pairs with energy asymmetry \E\ — E2\/{E\ + 
E2) < 0.7 is calculated and binned in pr- The 7r° yield 
in each pt bin is determined by integrating the back- 
ground subtracted two-photon invariant mass distribu- 
tion the background is determined from mixed events. 
The peak-to-background ratio increases from ~ 0.3 at 
pr=L25 GeV/c to more than 5 above 4.25 GeV/c. The 
raw 7T° spectra are corrected by Monte Carlo simulations 
for trigger efficiency, acceptance, and for 7r° reconstruc- 
tion efficiency including dead areas, effects of energy res- 
olution, photon identification cuts and peak extraction 
window. Finally, the yields are corrected to the center 
of the pt bin using the observed slope of the spectrum. 
Below pt— 5 GeV/c the yields are determined from the 
minimum bias data sample while above 5 GeV /c the pho- 
ton triggered sample is used. The main sources of sys- 
tematic errors arc listed in Table I. The final systematic 
errors on the spectra are 10 to 16%, increasing with pt- 

Charged particles are reconstructed using a drift cham- 
ber (DC) followed by two layers of multiwire proportional 
chambers with pad readout (PCI, PC3) In this 

analysis tracks were reconstructed over a restricted pseu- 
dorapidity range \r/\ < 0.18. Pattern recognition in the 
DC is based on a combinatorial Hough transform in the 
track bend plane, while the polar angle is determined by 
PCI and the location of the collision vertex along the 
beam direction |21| . The vertex was constrained to be 
within \z\ <18 cm for analysis of the charged tracks. The 
track reconstruction efficiency is approximately 98%, in- 
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TABLE II: Systematic errors in percent on the (h + + h )/2 
invariant yields. Error types are as in Table I. 





type 


pr <4 


Pt=4 


Pt=6 


PT=7 


PC3 match 


C 


2.4 


2.4 


2.4 


2.4 


geometr. acc. 


c 


2.9 


2.9 


2.9 


2.9 


Monte Carlo corr. 


c 


3.7 


3.7 


3.7 


3.7 


mom. resolution 


B 




<0.5 


0.6 


1.2 


mom. scale 


B 


<3.2 


3.2 


3.5 


3.7 


backgd. subtrac. 


B 


<0.5 


0.6 


3.8 


8.2 


7r oversubtrac. 


B 






1.8 


4.7 


Total error 




6.2 


6.2 


7.6 


11.5 


TABLE III: Systematic errors on 


the p+p reference spec- 


trum. 














type 


p T <2 


Pt=4 


Pt=G 


Pt = 10 


normalization 


C 


10 


10 


10 


10 


residual syst. 


B 


3 


5 


7 


9 


charged (ty°) only a 


B(C) 


8(5) 


7(5) 


9(5) 


20(5) 



"Errors affect only the charged (neutral) analysis, and include 
h/7ro, and fits to reference (point by point pp) errors. 



dependent of pt- Particle momenta are measured with 
a resolution Sp/p = 0.7% 9 l.l%p (GeV/c). The mo- 
mentum scale is known to 0.7%, from the reconstructed 
proton mass using the time-of- flight. A confirmation hit 
is required in PC3, located at a radius of 5 m, within 
a 2.5cr matching window to eliminate most albedo, con- 
versions, and decays. The remaining background above 
Pt=5 GeV/c is subtracted statistically using identified 
conversions and weak decays (primarily kaons) fisj . 

Corrections to the charged particle spectrum for geo- 
metrical acceptance, decays in flight, reconstruction ef- 
ficiency, and momentum resolution are determined us- 
ing a single-particle GEANT Monte Carlo simulation. 
All analysis steps, including the outer detector match- 
ing cuts are applied consistently in simulation and data. 
As the DC and PC occupancies in d+Au are low, 
no multiplicity-dependent occupancy corrections are re- 
quired. The yield was corrected to the center of the px 
bin. Table II lists the source and magnitude of each con- 
tribution to the systematic uncertainties on the charged 
particle spectra. 

The fully corrected px distributions of 7r°'s and (h + + 
h~)/2 are shown in Fig. Each panel shows the refer- 
ence spectrum from p+p collisions along with the par- 
ticle spectrum from d+Au. The reference for tt° is the 
7r spectrum measured in p+p by PHENIX shown 
on the right panel. For charged hadrons, the reference 
is obtained by scaling up the tt° spectrum by the pr- 
dependent h/V ratio observed at the ISR [22| at y/s = 23 
to 63 GeV and measured by PHENIX. The value of h/7r 
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FIG. 1: Midrapidity pr spectra for charged hadrons and 7r°. 
Total uncertainties are shown. The tt below 5 GeV/c arc 
from minimum bias triggered events, and above from photon 
triggered. Lines show fits to reference spectra from p+p col- 
lisions, and open points show the p+p it spectrum measured 
by PHENIX 



is constant at 1.6 ± 0.16 above 1.5 GeV/c, and decreases 
at lower px to a value of 1.25^q'25- Table III summarizes 
uncertainties on the reference spectra. 

A standard way to quantify nuclear medium effects on 
high pt production is by the nuclear modification factor, 
which we define for d+Au collisions as the ratio of invari- 
ant yield in d+Au to that of p+p, scaled by the number 
of binary collisions. 



R<ia(pt) 



{l/Nlf)d 2 N dA /dndp T 
(N coU )/a^ el <Pa pv /dr]dp7 



where (N co u) is the average number of inelastic nucleon- 
nucleon (NN) collisions per event in the minimum bias 
collisions, and (N co u) / cr™ el is the nuclear overlap func- 
tion (Ta(6)). Using a Glauber model and simulation 
of the BBC, (N co u) is 8.5 ± 0.4 in minimum bias d+Au. 

The ratio RdA is plotted separately for 7r° measured in 
the PbGl and in the PbSc calorimeters in Fig. [2] The two 
analyses are consistent within errors. The data are com- 
pared to the corresponding nuclear modification factor 
Raa obtained from central Au+Au collisions. The top 
panel of Fig. |3| shows RdA for inclusive charged particles 
(h + + h~)/2, again compared with Raa observed in cen- 
tral Au+Au collisions, while the lower panel compares 
(h+ + h-)/2 with tt . In both Figs. Eland the uncer- 
tainties are plotted as follows: error bars represent the 
quadrature sum of statistical errors and those systematic 
errors which vary point-to-point in pt; systematic errors 
on the absolute yield and the systematic errors which are 
correlated point-to-point are shown as bands. All proce- 
dures for estimating the various systematic errors have 
been adjusted to provide estimates corresponding to la 
error values. 
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FIG. 2: Nuclear modification factor R^a for it in the PbGl 
and PbSc calorimeters in minimum bias d+Au. The bands 
around the data points show systematic errors which can vary 
with pt, while the shaded band around unity indicates the 
normalization uncertainty. The nuclear modification factor 
Raa in 10% most central Au+Au collisions is also shown. 

The data clearly indicate that there is no suppression 
of high pt particles in d+Au collisions. We do, how- 
ever, observe an enhancement in inclusive charged parti- 
cle production at px > 2 GeV/c. A similar enhancement 
was observed in p+A fixed-target experiments |24| and 
is generally referred to as the "Cronin effect" . To facili- 
tate comparison of the Cronin effect in inclusive charged 
particles and 7r°, the lower part of Fig. [3] shows all sys- 
tematic uncertainties common to both analyses in the 
bar on the left. It should be noted that this uncertainty 
must be added in quadrature with the bands shown for 
each curve to obtain the la allowed range of RdA from 
the data. The ir° data suggest a smaller enhancement 
for pions than for inclusive charged particles at pr = 2- 
4 GeV/c. We note that the charged spectrum includes 
baryons and antibaryons, which may have a different nu- 
clear enhancement than the mesons |24| . 

The various models of the suppression observed in 
Au+Au predict a different dependence on N co u in 
d+Au 0, |2ji |2(|. Therefore, a second data sample 
was selected by requiring observation of a neutron in the 
Zero-Degree Calorimeter on the deuteron-going side of 
PHENIX. This, together with the requirement of parti- 
cles entering both Beam-Beam Counters, selects a class 
of events in which only the proton from the deuteron in- 
teracts with the Au nucleus. The mean number of binary 
collisions for this sample is calculated with the Glauber 
model to be 3.6 ± 0.3. Particle yields in this sample 
have a < 5% uncertainty beyond that of the minimum 
bias sample, arising from trigger bias. 

Fig.0|shows the ratios of RdA in minimum bias d+Au 
to R p a in the neutron tagged sample, for both (h + + 
h~)/2 and ir° . Systematic uncertainties on the spectra 
cancel in the ratio; the band around unity shows the 



05 2 F 




\ r n 

0123456789 10 

p T (GeV/c) 



FIG. 3: Top: Nuclear modification factor .RdA for (h + +h~) /2 
in minimum bias d+Au compared to -Raa in the 10% most 
central Au+Au collisions. Inner bands show systematic errors 
which can vary with pr, and outer bands include also the 
normalization uncertainty. Bottom: Comparison of -RdA for 
(h + +h~)/2 and the average of the tt measurements in d+Au. 
The bar at the left indicates the systematic uncertainty in 
common for the charged and ir° measurements. 



uncertainty on the ratio of the number of binary collisions 
in the two samples. Average values of N co ii are 3.6 per 
participating proton in the neutron tagged sample and 
8.5 for 1.7 participating nucleons from the deuteron in 
minimum bias d+Au. Given the systematic uncertainties 
on N co u , we cannot exclude a small centrality dependence 
for pt > 1 GeV/c. It should be noted that the figure also 
indicates that d+Au collisions provide a good measure of 
the physics of p+Au. 

The observation of an enhancement of high-pr hadron 
production in both the minimum bias d+Au and the neu- 
tron tagged sample of p+Au collisions indicates that the 
suppression in central Au+Au collisions is not an initial 
state effect. Nor does it arise from modification of par- 
ton structure functions in nuclei. The data suggest, in- 
stead, that the suppression of high pt hadrons in Au+Au 
is more likely a final state effect of the produced dense 
medium. 
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FIG. 4: Ratio of _RdA in minimum bias d+Au and R p a from 
neutron tagged d+Au collisions for (h + + h~)/2 and 7r°. The 
band shows the uncertainty on the number of binary colli- 
sions; all other systematic errors cancel. 
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